SEMI-CLASSICAL LIMIT OF SCHRODINGER POISSON 
EQUATIONS IN SPACE DIMENSION n > 3 



THOMAS ALAZARD AND REMI CARLES 



Abstract. We prove the existence of solutions to the Schrodinger- 
Poisson system on a time interval independent of the Planck constant, 
when the doping profile does not necessarily decrease at infinity, in the 
presence of a subquadratic external potential. The lack of integrability 
of the doping profile is resolved by working in Zhidkov spaces, in space 
dimension at least three. We infer that the main quadratic quantities 
(position density and modified momentum density) converge strongly as 
the Planck constant goes to zero. When the doping profile is integrable, 
we prove pointwise convergence. 



1. Introduction 

We consider the semi-classical limit e — > of the Schrodinger-Poisson 
system: 

£ 2 

(1.1) ied t u e + — Au £ = V ext u £ + V £ u £ , (t, x) e E x R n , 

(1.2) AV£ = q{\u e \ 2 -c), 

(1.3) uf t=0 = alixY^'^ 

where V ex t = V ex t(t,x) is an external potential (harmonic potential for in- 
stance), c = c(x) is a doping profile (or impurity, background ions), and 
q £ M represents an electric charge; V ex t, c and q are data of the problem 
(see e.g. |2fl|). We consider the case where the space dimension is n ^ 3. 
This is due to a lack of control of low frequencies for the Poisson equation 
(fj~2|) when n < 2. 

The conditions we impose to solve the Poisson equation (|1.2j) will be given 
according to the different cases we consider. 

The doping profile c is supposed to be bounded, and does not neces- 
sarily goes to zero at infinity (see Assumption ^ or Assumption [21 below). 
Suppose for instance that c = 1. Then ()l.lj) - (|1.2j) is reminiscent of the 
Gross-Pitaevskii equation (see e.g. ^BlEI] an d references therein): 



(1.4) ied t u £ + — Au £ = (\u £ \ 2 - 1) u £ . 
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For this equation, the Hamiltonian structure yields, at least formally: 
|(|| £ V^(t)||i 2 + |||n £ (t)| 2 -l||^) =0. 

A natural space to study the Cauchy problem associated to (|1.4[) is therefore 
the energy space 

E = {ue Hl c (R n ) ; Vu G L 2 (R"), \u\ 2 - 1 G L 2 {R n )}. 

For this quantity to be well defined, one cannot assume that u £ is in L 2 {R n ); 
morally, the modulus of u e goes to one at infinity. To study solutions 
which are bounded, but not in L 2 (R n ), P. E. Zhidkov introduced in the 
one-dimensional case in [25] (see also |26j ) : 

X s {R n ) = {u G L°°(R n ) ; Vu G H s ^ 1 (R n )}, s > n/2. 

The study of these spaces was generalized in the multidimensional case by 
C. Gallo |IIJ. They make it possible to consider solutions to (|1.4|) whose 
modulus has a non-zero limit as |x| — ► oo, but not necessarily satisfying 
\u £ (t,-)\ 2 - 1 G L 2 (R n ). 

Recently, P. Gerard JB] solved the Cauchy problem for the Gross-Pitaevskii 
equation in the more natural space E, in space dimensions two and three. 
The main novelty consists in working with distances instead of norms, in 
order to apply a fixed point argument in E. In particular, the constraint 
\u e (t, -)| 2 - 1 G L 2 (R n ) is satisfied (and propagated). 

We have to face a similar issue, when solving the Poisson equation. Mim- 
icking the approach of ^2E3) ^ * s na tural to work with the property: 

\u £ (t,.)\ 2 -c(-)€L 2 (R n ). 

We shall always assume that this holds at time t = 0. We prove that this 
property holds on [0, T] for some T > independent of e, provided that we 
consider an external potential whose unbounded part is linear in x. However, 
our analysis shows that in the presence of a quadratic external potential, this 
property is not relevant off t = (see Section [SJ. 

Note that we make no assumption on the sign of q (which models the 
charge of the element considered in a semiconductor device). This is in 
sharp contrast with the mathematical analysis of the semi-classical limit 
of the nonlinear Schrodinger equation. When the Poisson term Vp(t,x)u £ 

is replaced with the nonlinear term f(\u £ \ 2 )u £ , E. Grenier ^H] proposed a 
strategy to obtain a phase/amplitude representation of the solution u £ . This 
leads to study a quasi-linear system whose principal part writes: 

D f , :=d f 2 -div(/'(H 2 )V-). 

Hence, to prove that the Cauchy problem is well-posed, one has to assume 
that the nonlinearity is defocusing and cubic at the origin (/' > 0), except for 
analytic initial data ^2], for which one can solve elliptic evolution equations. 
Here, we are not restricted to the case when q > 0. As will be clear below, 
the reason is that the quasi-linear operator Hf is replaced with the semi- 
linear operator d 2 — gA -1 V((|u £ | — 1) div •). 
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Notation. Recall that for s > n/2, Zhidkov spaces are defined by 1 : 
X s (R n ) = {ue L°°(R n ) ; Vu G F^R™)}- 

We denote 

We write H s = H s (R n ), X s = X s (R n ), H°° := n seN # s , X°° := n seN X s . 
We do not use specific notations for vector-valued functions: for instance, 
we write abusively V 2 / G H°° when d? k f G H°° for every 1 ^ j, k ^ n. 

Remark 1.1. Zhidkov spaces contain all the functions of the form 

7 + v, with 7 = Const. G C and v G H s (R n ). 
The converse is not true, as shown by the following example: 

u ( x ) = -. X ) n ' X = (xi,X2,X 3 ) eR 3 . 

1 + \x\ z 

On the other hand, if n ^ 3 and u G X s for some s > n/2, then there exists 
7 G C such that u — 7 G L n ~ 2 (R n ) (see Lemma 12 . 1 1 b elow) . 

In this paper, we consider the system (jl.lJI - Ql.3JI in three cases: 

• The external potential and the initial phase are sub- linear in x, and 
the mobility c is in Zhidkov spaces (Part ^) . 

• The external potential and the initial phase are sub-quadratic in x, 
and c is a short range perturbation of a non-zero constant (Part [2]). 

• The mobility is integrable, and the external potential and the initial 
phase are sub-quadratic in x (Part EJ) . 

In the first two cases, we construct a solution to Q1.1J1 - Q1.3J) in Zhidkov 
spaces, and describe the asymptotic behavior of the main quadratic observ- 
ables as e — > 0. In the last case, we construct a solution in Sobolev spaces, 
and give pointwise asymptotics of the solution as e — > 0. 

In this introduction, we describe more precisely the results corresponding 
to the first case. We emphasize the fact that if we simply assume V eK t G 
C(R;H°°) and $0 £ H°°, then our analysis becomes much simpler. The 
unboundedness of V ex t and <3?o require some geometrical description that 
complicates the technical approach. Yet, this makes our assumptions more 
physically relevant (see e.g. ^1] and references therein). 

Assumption 1. Recall that n 3. 

• External potential: V ext G C°°(R x R n ) writes 

V cxt {t,x) = E(t)-x + V pevt {t,x), with E G C°°{R) and W port G C(R; H°°). 

• Doping profile: c G X°° . 

• Initial amplitude: a^{x) = do(x) + r £ (x), where a$ G X°° is such that 
|d | 2 - c G L 2 (R n ), and r £ G H°° , with 

IKHffs — >0, Vs 5s 0. 

• Initial phase: we have $0 £ C°°(lR n ) with 

$o(x) = a ■ x + (f)o(x), with a G R n and V<fio G H°°. 

^Fot general s > 0, another definition is used, see 
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Lemma 1.2. Under the Assumption^ there exists a unique solution e ik 6 
C°°(]R x R n ) to: 

(1.5) d^eik + ^|V(/> eik | 2 + £(£) -x = ; </> ei k(0, x) = a ■ x + O ■ 
This solution is given explicitly by (f)^{t,x) = ait) ■ x + (3{t), where: 

a(t) =a - [ E(r)dT ; (3{t) =(3 Q - \ [ air) 2 dr. 
Jo 2 J 

We skip the proof of this lemma; a more general result is proved in Sec- 
tion El We will see that if V ex t and/or $o have a quadratic dependence on 
x, then we have to consider an eikonal phase ^ e ik which is quadratic in x. 

Theorem 1.3. Let Assumption^ be satisfied. There exists T > indepen- 
dent ofe G]0, 1] and a solution u £ G L°°([0,T] x W 1 ) to (fTT]l - (fOj) . with 

VV £ {t, x) -> as \x\ -> oo, F p e (i, 0) = 0, 

and such that \u £ \ 2 — c G L°°([0, T]; L 2 ). Moreover, one can write u £ = 
a £gi(</) e i k +0 e )/e^ where: 

• a £ G C°°([0, T] xR n )nC([0,r];X°°), and \a £ \ 2 - c G C([0, T]; L 2 ). 

• 4> £ £ C°°([0,T] x M n ) andV(p £ G C([0,T];X°°). 

• We /wive i/ie following uniform estimate: for every s > n/2, there 
exists M s independent of e G]0, 1] such that 

\\ a£ \\L°e(0,T;X s ) + IlKI 2 _ C \\l°°(0,T;L 2 ) + \\l°°(Q,T;X s ) ^ Mj- 

Remark 1.4. We could not prove a uniqueness result for u £ . 

Remark 1.5. The above conditions to solve the Poisson equation are similar 
to those given in |24j . We explain at the end of Section 13.31 why in our 
framework, we cannot impose V £ (t,x) — > as \x\ —* oo (as in [21 HS1 f° r 
instance) . 

Besides the uniform bounds, even the existence of such a solution u £ 
is new. First, the presence of the external potential seems to have never 
been studied rigorously before. As we already mentioned, this makes the 
proof more technically involved. Next, in most of the previous studies, u £ 
is supposed to be in L 2 : see e.g. |^. In (23], the author considers 
the case c G L 1 n H s . As we will see in Section |S1 this case makes the 
analysis easier, and also makes it possible to have u £ G L 2 . The main 
difficulty in the analysis lies in the fact that when c G" L 1 (M. n ), the condition 
\u £ \ 2 — c G L 2 (R") is somehow "more nonlinear", as in |13| . 

The general idea to prove Theorem 11.31 consists in adapting the idea of 
|15j : with techniques from the hyperbolic theory, we construct a solution to 



1 

T 

1 £ At r £ _ ■£ 

-a A<£ — ^2"" ' "|t=0 — "0- 
AV £ = q(\a £ \ 2 - c) ; VV £ (t,x) — > 0, V £ (t,0) = 0. 



d t <S> £ + ^\V$ £ \ 2 + V ext + V £ = ; $f t=0 = $o- 
(1.6) d t a £ + V$ e • Va £ + -a £ A<S> £ = i-Aa £ ; al_ n = af 
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Following |5j, we write $ e = (f) e ^. + 4> £ : with the unknown (a £ ,(f) £ ), 1)1.6(1 
becomes (we keep the term A0 e ik which is zero here, for future references): 

d t <p £ + V^ik • V<f + ^|V<f| 2 + V pert + V p e = 0; 0f t=o = O . 

(1.7) d t a £ + V (<f + 0eik) • Va £ + ^a £ A (<f + eik ) = i^Aa £ ; af t=0 = ag. 
AV p e = g(|a e | 2 - c) ; VV £ (t,x) -— 0, V p £ (i,0) = 0. 

Proving the existence and uniqueness of solution to (|1.7|) as we do in the 
proof of Theorem ll.3l is one of the main results of this paper. Because of the 
difficulties pointed out above, and the fact that one can easily be mistaken 
by using the usual approach, we give full details for the construction of the 
solution to p. 7)1 . Passing formally to the limit, it is natural to consider: 

d t (f> + V0 cik • V<P + -|V0| 2 + Vpert + V p = ; 0| t=o = O . 

I 

(1-8) <9 t a + V (0 + cik ) • Va + ~aA (0 + eik ) = ; a\ t=Q = a . 

AV P = g(\a\ 2 - c) ; VV p (t, x) — > 0, V p (t, 0) = 0. 

\x\^oo 

Notation. The symbol < stands for ^ up to a positive, multiplicative con- 
stant which depends only on parameters that are considered fixed. 
We shall also denote L™Y for L°°([0,T];Y). 

Theorem 1.6. Under Assumption^ there exists a smooth solution (a, 0) 
of (JUJ) such that a, V0 G C([0, T], X°°), |a| 2 - c e C([0,T],L 2 ), and 



a||L2?tfs + ||V(0 e - 0)||z°fA- — >0, Vs > n/2. 

- 1 J e^0 



7n particular: 



\u £ \ 2 — Hal 2 inL^H 3 , and 

elm(WVu £ ) — >|a| 2 V(0 eik + (/>) inL^X 3 , Vs > n/2. 

Recall that in general, none of the terms a or a e is in L 2 (M. n ). Though, 
the difference a 6 — a is in L 2 (M n ), and asymptotically small as e — > 0. Note 
that (p, v) := (|a| 2 , V(0 + </> c ik)) solves the Euler-Poisson system: 

<9 4 p + V • (pv) = 0, 
&v + v • Vv + VV ext + W D = 0, 



(1.9) 



'est "I" V V p 

AV P = q(p-c) ; W p (t, x) — > 0, F p (t, 0) = 0. 

|at|— »oo 



The existence of solutions to p. 9(1 under Assumption ^ is new. 

This paper borrows several ideas from [5J, and [T^l- As we have al- 
ready mentioned, an important difference with |15j is that the underlying 
wave equation associated to p.6|) is semi-linear, and not quasi-linear. The 
reduction to p.7[) is similar to the approach in [5J- Several important dif- 
ferences should be pointed out. First, we work in Zhidkov spaces instead of 
Sobolev spaces, an aspect which requires some extra care. Integrating the 
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Poisson equation, especially when we have AV^f G L 2 (M. n ) and not neces- 
sarily AVp G is also a new problem. Finally, the propagation of 
the initial assumption |q,q| 2 — c G L 2 (M n ) turns out to be different from the 
phenomenon studied in JH] • As we shall see in Section 03 the presence of 
quadratic "geometric" quantities (such as an external harmonic potential) 
requires a highly non-trivial adaptation of the approach in jSJ. 

The rest of this paper is organized as follows. In Section [2j we collect 
various technical estimates, in order not to interrupt the proofs later on. In 
Section EJ we prove Theorem 11.31 Theorem 11.61 is proved in Section ^ In 
Part [21 (Sections EH3)) we consider the case when c — 1 G L 1 n H°°, and the 
external potential and the initial phase contain quadratic terms. In Part 01 
(Section^, we assume c G L 1 (M™), and prove a refined convergence result. 

Remark 1.7. Before leaving this introduction, let us explain why we concen- 
trated on the whole space problem. Indeed, some problems require consid- 
ering the periodic case (see [Q and the references therein), where the space 
variable belongs to the torus T n . As a matter of fact, the periodic case is 
easier. This follows from two observations: first, the computations below 
apply mutatis mutandis in the periodic setting; and second, for all a G M, 
the operator A _1 V is well-defined in H a (T n ). 

2. Estimates in Lebesgue, Sobolev and Zhidkov spaces 

This section serves as the requested background for what follows. The 
proofs of easy or classical results are left out. We first recall a conse- 
quence of the Hardy-Littlewood-Sobolev inequality, which can be found in 
[TCI Th. 4.5.9] or 13. Lemma 7]: 

Lemma 2.1. // cp G £>'(IR n ) is such that V<p G L p (W n ) for some p G]l,n[, 
then there exists a constant 7 such that ip— 7 G L q (M. n ), with 1/p = 1/q+l/n. 

This shows that under Assumption ^ the doping profile is of the form 

c = 7 + c, where 7 is a constant, and c G L™- 2 (R™), Vc G H°°. 
Define the Fourier transform as 

Lemma 2.2. Let n ^ 3. For every s > n/2, there exists C s such that 

(2.1) \\y\\ L °° ^ CjV^lff.-!, V93 G H s (R n ). 

Remark. In space dimension n ^ 2, low frequencies rule out the above 
inequalities. For instance, in space dimension n = 1, the function 

f( x ) = f = ar g sinh(x) 

Jo vi + y 

is not in L°°(M), but its derivative is in H°°. In space dimension n = 2, 
consider the function 

f{x 1 ,x 2 ) = log |log(xf + . 
One can check that V/ G H°°, while clearly, / L°°(R 2 ). 
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Warning (Homogeneous Sobolev spaces). It may be tempting to restate 
Lemma l'2.2l in terms of homogeneous Sobolev spaces. Recall that, for s > 
0, the homogeneous Sobolev space H s is defined as the completion of the 
Schwartz space 5(IR n ) for the norm 

More precisely, one might want to replace the right-hand side of (|2.1j) with 
IMIiji + IMIijs an d consider tp G H 1 n H s only. This is extremely delicate, 
since H s is not a Hilbert space when s ^ n/2. 

Lemma 2.3. Let n ^ 3, q ^ 2 and s > n/2 — 1. There exists C = C(n, q, s) 
such that for all tp G L q (R n ) with Vip£ H s (R n ), 

\\ip\\ L oo ^C(\\ip\\ Lq + \\Vp\\ H s) . 
Proof. The usual Sobolev embedding yields, for any a > n/q, 

IMIl°° ^ + II|v| <t v'||l9. 

On the other hand, for k = n(l/2 — 1/q), 

iim'viil. < mvr^ii^ < iiv^h^-!, 

provided that a ^ 1. If s > n/2 — 1, a given by s = n(l/2 — 1/q) + a — 1 
is such that a > n/q and <r ^ 1. The above two estimates then yield the 
lemma. □ 

Lemma 2.4. Let n ^ 3. For every s ^ 0, VA" 1 maps L 1 (R n ) n H s (R n ) to 
H s+1 (R n ): there exists C s such that 

||VA -1 y>[| H .+i C s (l^llii + , G L 1 (R n ) n # s (K n ). 

The following variant of the classical Kato- Ponce estimates can be found 
in |17| Theorem 5]: 

Lemma 2.5. Let n ^ 1 and s > n/2 + 1. Denote A = (L — A) 1 / 2 . T/iere 
exists a constant C s such that, for all f G X s+1 (M n ) and all u G H s ~ 1 (W l ), 



(2.2) ||/A s n - A s (/n)|| i2 < C s (||V/|| xoo |M| HS -a + HV/H^ ||«|| roo ) . 

Lemma 2.6. Lei s > n/2. The Sobolev space H s (R n ) and the Zhidkov 
space X s (R n ) are algebras: there exists a constant C s such that, for all 
u,v G H s (R n ) and a,b £ X s (R n ), 

\\UV\\jj s 

u \\h s \\ v \\h s i ll^llx 8 ^ C s || a |lx s ll^llx s 1 

There exists C s such that for all v G H s (R n ) and a G X s (R n ), 

\\av\\ Hs + ||au|| xs < C s \\v\\ H s \\a\\ X s ■ 

There exists C s such that for all a G X s (R n ) and b G X s+1 (R n ), 

\\a7b\\ H s < C s \\a\\ X s \\b\\ X s+i ■ 

In order to use Arzela-Ascoli's theorem, we will invoke: 

Lemma 2.7. Let a > n/2 and ((fj)j^ be a bounded sequence in X a (R n ). 
For all a' < a, there exists a subsequence which converges in (W 1 ). 
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Proof. This follows from the fact that, for all test function \ G Co°(R n ), 
(x^j)jeN is a bounded sequence in H a (W l ). □ 

Remark. It might seem more natural to state a precompactness result in 
Z£' C (R") := [p G L£" C (R") ; G ^(R")} ■ 

Actually, one can check that for a' > n/2, X^ c (R n ) = H^ c (R n ). 

Lemma 2.8. Let n ^ 3 and s > n/2. 

• For all p > ^zo, there exists C = C(s,p,n) such that: 

IIA^V/Ulp <C\\f\\ H ; VfeH s . 

• There exists C = C(s,n) such that: 

Proof. Essentially, we use the property / G L 2 for low frequencies, and 
/ G L 1 for high frequencies (J-(H S ) C L 1 if s > n/2). For p > 2n/(n — 2), 



< 



LP 



L^F(]C|<1) 



L 2 



L°°(|£|^l) 



The norms involving |£| 1 are finite since p > 2n/(n — 2). For s > n/2, 



< 



(0 s f 



L 2 



The first point follows from the Hausdorff- Young inequality: 



2/ 



Lp' 



The second point is straightforward, with p' = 1. 



□ 



Part 1. Sublinear eikonal phase 

3. Proof of Theorem 11.31 

Our first task is to construct a solution to Q1.7JI . As explained in the in- 
troduction, it is convenient to introduce the "velocity" v £ = V(p £ . Denoting 
v eik = V0 e ik, and recalling that is a function of time only, we infer from 
(|1.7|) that (a £ ,v £ ) has to solve: 

( 8 t v £ + (v eik + w £ ) • W + w pcrt + vv £ = 0, 



(3.1) 



together with 
(3.2) VV £ (t,x 



d t a £ + (v cik + v £ ) ■ Va £ + ^a e V • v £ = i^Aa £ , 
I AF p £ = g(|a £ | 2 -c), 

; V £ (t,0) = ; v? t=0 = V0 O ; af t=0 = a§. 



|a;|— >oo 



In the context of Assumption ^ we show that the solutions of (|3.1j) - (|3.2j) 
exist and are uniformly bounded for a time interval independent of e. 
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Proposition 3.1. Let Assumption^ be satisfied. Let s > n/2. For all 
M > Mq > ; there exists T > such that, if for all e G [0, 1], 

(3.3) ||V<fo||# s +2 + ||K| 2 - c\\ L2 + \\a £ \\ X s+i < Mq, 

then the Cauchy problem (|3.1|) - (|3.2|) has a unique classical solution (v £ ,a £ ) 
in C°°([0,T] x W 1 ) such that 

(3-4) \\v £ \\ L¥X s+2 + |||a £ | 2 - c|| L??L2 + \\a e \\ L¥XS +i ^ M. 

As suggested by the above statement, we construct W £ (only the gradi- 
ent of V £ is present in (|3,1|) ), and the condition V £ (t, 0) = is given only to 
insure uniqueness for V £ (even though it is not stated in the above result). 
Therefore, we shall neglect this condition for a while. 

3.1. Regularized equations. Let j be a C°° function of £ G R n , with 

0<J<1, j(0 = 1 for |£K 1, XO = for |£| ^ 2, j(C)=j(-0- 
Set jh(£) := j(h£), for /i > and £ G M n ; jh is supported in the ball of radius 
2/h about the origin. Define Jh as the Fourier multiplier with symbol jh- 

Jh ■= j(hD x ). 

Also, for our purpose it is interesting to introduce a family of operators 
that cut the low frequency component of a function. Indeed, the Poisson 
term W p e = qA^V^l 2 - c), is not well defined in general. We replace 
the operator (/A _1 V by a family of operators i?hV well defined on Sobolev 
spaces and prove that, in the end, there is no need to estimate the low 
frequency component of W £ . To do that, we set 

Gh = I - J\/hi 

that is, Gh is the Fourier multiplier with symbol 1— Ji/h, which is supported 
in {|£| ^ h}. Consequently, the operator 

R h := qA^Gh, 

is bounded in all Sobolev spaces (with operator norm going to +oo when 
h tends to 0). More precisely, there exists a constant C such that, for all 

° > °> 

\A G}i 1 1 > ^ja+2 ^ Ch . 

Consider the following approximation of (|3.1|) : 

( d t v e h + J h {(v cik + v £ h ) ■ VJ h v £ h ) + W pcrt = -R h V{\a £ h \ 2 - c), 

(3-5) < 1 £ 

[ d t a% + J h ((v eik + v £ h ) ■ VJ h a £ h ) + -a £ h V ■ v £ h = i-AJla £ h . 

We keep the same initial data: 

( 3 - 6 ) v l\t=o = V( ^o ; a%\t=o = a o- 

Note that Assumption ^ implies that v^ t=0 is in and is independent of 
e G [0, 1] and h G]0, 1], while a|u =0 is in X°°, and uniformly bounded in X s 
for any s > n/2, for e G [0, 1] and h G]0, 1]. 

The point is that the regularized equations (|3.5|) ~ (|3.6|) nave been chosen 
so that the Cauchy problem can be solved as in the standard framework of 
Sobolev spaces: 
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Lemma 3.2. Let s > n/2. For all e G [0, 1] and all h E]0, 1] there exists 
> suc/i i/iai the Cauchy problem Q3.5j) - Q3.fij) /ias a unique solution 
(v £ h ,a £ h ) E C 1 ([0,2£];fr a+2 (R n ) x X s+1 (M n )). 

Proof. The proof is based on the usual theorem for ordinary differential 
equations. Set u £ h = (v^a^) and we rewrite Q3.5j) under the form 

d t u% = F x {e, h, u%) + F 2 (t)u £ h + F 3 (t, x), 

where Fi(e, h, u) is at most quadratic in u, and we have used the property 
that f e ;k is a function of time only. We have to verify that the functions F 
are smooth. This follows from Lemmas 12.21 and 12.61 and the fact that the 
operators Rh and AJh are of order —2 and respectively: 

\\J h (v £ h ■ VJ h V £ h )\\ HS +2 < IKH.H-,+2 \\VJ h V £ h \\ H s+2 < hT X , 

\\Jhi v h ■ VJha 6 h)\\x s + 1 ~ \\^Jha £ h\\xs+ 1 ~ h l a^+i > 

\\ a fy ' v h\\x s + 1 ~ ll a /illx s + 1 ' 

\\RhV\a £ h \ 2 \\ HS+2 </i- 2 ||VK| 2 || ff . <h- 2 \\ai\\ 2 X3+1 , 
\\ AJ h a h\\ x *+i = \\ J hAJ h a e h \\ X s+i ^ \\AJ h a%\\ HS+ i < h~ 2 \\a £ h \\ xs+1 . 

□ 

3.2. Uniform bounds. To prove Proposition 13. 11 the analysis of (|3.5|) con- 
tains at least two parts: first, an existence and uniform boundedness result 
for a time independent of the small parameters e and h; and second, a con- 
vergence result when h — ► 0. Here, we prove that the solutions (v^, afj exist 
and they are uniformly bounded for a time independent of the parameters 
e and h. Below, T^* denotes the lifespan, that is the supremum of all the 
positive times such that the Cauchy problem for H3.5jl - Q3.fijl has a unique 
solution in C^fO, T^); H s+2 (R n ) x X s+1 (M n )). 

Proposition 3.3. Let s > n/2. There exists a continuous function g: R* + — > 
R* + such that, for all e E [0, 1] and all h e]0, 1], the norm M £ h : [0, Ty*\ -> R* + 
defined by 

M h( T ) : = \\ a h\\L^X^ + \\\ a h\ 2 ~ c ILf?L 2 + l|V^|| L oc> Hs+ i, 

satisfies the estimate: M £ (T) ^ M £ (0)e T 9( M h( T )\ VT G [0,T £ *]. 

Proof. Before we proceed, two comments are in order. Firstly, the functions 
(v £ h ,a £ h ) are smooth (C 1 in time with values in Sobolev/Zhidkov spaces), so 
that it is easily verified that all the following computations are meaningful. 
Secondly, it is useful to note that, in view of Lemma l2~2l it suffices to prove 
that 

(3.7) m £ h (T) ^ m£(0)e T *W)), VT G [0,Tf ], 

where 

= ||Va|||i,|?i^ + ||K| 2 - c|| L§?L2 + ||Vt^|| Z oo H *4-i, 

and 

Mh(T) := Mfi(T) + ||^IIz,°o([o,t]xR«) • 
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Indeed, Lemma 12.21 provides us with a constant C s such that ^ C s Mf, 
and we have: 

Lemma 3.4. Let s > n/2 and c G X°° . There exists a constant K such 
that, for allT > and if £ X°° , 

IMli°°(R«) ^ K \\\f\ 2 ~ c\\ L 2 + K\\V(p\\ 2 H s +K\\c\\ X s+i. 

Proof. By Lemma 12.31 we have: 

||M 2 - c\\ LX < ||M 2 - c|| L2 + ||v(M 2 - c) H^. 

Since s > n/2, 

||V(M 2 - c)\\ H . < M Loa \\V<p\\ H . + \\V<p\\ 2 Hs + \\Vc\\ Hs . 
Triangle inequality yields 

\\<p\\l<*> ~ \\\f\ 2 ~ c IIl2 + ll^llioo ||V<p|| Hs + ||V^||^ S + ||c|| xs+1 , 
hence, the desired result follows by Young's inequality. □ 

With these preliminaries established, to prove (|3.7|) . we begin by estimat- 
ing the L 2 norm of \a 6 h \ 2 — c. To do that, we start from 

^IIKI 2 - 4 2 L2 < 2|||4| 2 - c|| L2 ||^(|4| 2 - c)\\ L2 . 

The second factor in the right hand side is estimated by 
||<MK! 2 -c)|| L2 < 2||o^|| LOO \\d t a £ h \\ L2 . 
Directly from the equations, we find that for bounded times, 

\\dta £ h \\ L2 < (1 + H\\ LOO ) \\Val\\ L2 + \\a%\\ Loa \\Vv e h \\ L2 + \\Aa £ h \\ L2 . 
Consequently, we obtain 

(3-8) j t \\\al\ 2 -c\\ 2 L2 <MUl + M £ h ) 2 . 

We now turn to the estimate of the H s norm of Vaf . Set Q := A S V, where 
A = (J - A) 1 / 2 . Since [V,Q] = = [J h ,Q], by commuting Q with the 
equation for a £ h , we find: 



with 



d t Qa% + J h ({v cik + v £ h ) ■ VJ h Qa £ h ) - i~AJ£Qa% = f £ h , 



f £ h := J h ([v%, Q] ■ VJ h a £ h ) - l -Q{a £ h V ■ v £ h ) 



Notice that is self-adjoint. We use the following convention for the scalar 
product in L 2 : 



(f,9)-= / f{x)g{x)dx. 
We have, since Vv e ik = 0: 

Re( iAJ^Qal , Qa £ h ) = Re( iAJ h Qa £ h , J h Qa% ) = 0, 
2( M(veik + v%) ■ VJ h Qa £ h ) ,Qa £ h ) = 2( {v cik + v%) ■ VJ h Qa £ h , J h Qa £ h ] 

= -({V -v £ h )J h Qa £ h , J h Qa%). 
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Therefore, 



f t \\Qa%f L2 =2Re(d t Qal,Qat) 

= Re( (V • v £ h )J h Qa £ h , J h Qa% ) + 2 Re( f h , Qa\ > 

^ IIV^H^ \\JhQa e h \\ 2 L2 + 2 \\m\ L 2 \\Qa%\\# . 

We now have to estimate the L 2 norm of The first term is estimated by 
way of the commutator estimate (12.21) and the Sobolev embedding: 

\\J h ([v E h ,Q]-VJ h al)\\ L2 < \M,Q]-VJ h al\\ L2 

< (l|Vt£|| LO o + llV^H^+x) \\VJ h a%\\ Hs 

< llV^H^+x \\Val\\ Hs . 
To estimate the last term, we use Lemma 12.61 to obtain 

\\Q[aiy-v E h )\\ L2 < \\a%V ■ vl\\ He+1 < \\a%\\ xs+1 ||V^|| HS+1 . 
We infer that 

11/hll.L 2 ~ II^aIIb-'+i ll a ft,llx s + 1 • 

Therefore, we end up with 

(3-9) l l|Va ^^ s ~ \\VvUhs+i\K\\xs+i- 

The technique for estimating Vu? in H s+1 is similar. Indeed, the analysis 
establishing the previous estimate also yields 

|||V^||| S+1 < (i + WVvMh^) \\VvUW 



+ || W pert ||^ + i + HV^Vdall 2 - c) \\ 2 H 



s + l • 



Since VR^V is uniformly bounded from H s+1 to itself, we obtain 

|||V^||^ +1 < (i + ||v^|| H3+1 ) \\Vvi\\W + 1 + ||I4| 2 - cf HS+1 . 

Next, noting that 

\W a h\ 2 ~ c |li^+i ~ \W a h\ 2 ~ C || L 2 + ||Vc||h« + l|V|a|| 2 ||^ 

~ IlKI 2 - c ||l2 + 1 + IKHi- + ll v °lll^ 
^ ||l a ll 2 ~~ c IL 2 + 1 + II ci /iIIa' s + 1 ' 

we conclude that 

(3-10) j t \Wv £ h \\ 2 HS+1 < C{M £ h ) || Vv%\\ 2 HS+1 +C(M £ h ). 

Summing over ()3.8|) . ()3.9j) and (|3.1U|) . Gronwall lemma yields the uniform 
estimate (|3.7j) . □ 

Lemma 13.21 and Proposition 13.31 yield the following result: 

Corollary 3.5. Let Assumption^ be satisfied, and let s > n/2. For all 
M > Mq > 0, there exists T > such that, if for all e € [0, 1], 

||V(/>o||#s + 2( R n) + |||ao| 2 ~~ C (')||i2( R n) + 1 1 °0 1 1 A^ 1 (R™ ) ^ ) 
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then the Cauchy problem ()3.5j) - (j3.6j) has a unique classical solution (y?, af) G 
C 1 ([0,T];fl' s+2 (]R n ) x A s+1 (M n )) satisfying 

\\Vv £ h \\ L cv H s+i + \\\a e h \ 2 - c|| L oo L 2 + ll a /illL5?x s + 1 ^ M - 

Remark 3.6. Refining the above computations thanks to Moser's calculus 
and tame estimates, we can see that the above existence time T can be taken 
independent of s > re/2 (see e.g. ^3 Section 2.2] or [^H Section 16.1]). 
This explains why we did not emphasize its dependence upon s, and why 
we consider different values for s below, without changing the notation T. 

3.3. Convergence of the scheme. We first claim that d^vf and dta £ h are 
bounded in C([0,T]; X s " 1 ), uniformly for h G]0, 1]. To see this, by us- 
ing Lemma 12.21 ()3.5|) and Corollary 13.51 the point is to verify that the 
term RhW (|a|| 2 — c), in the equation for dtv £ h , is uniformly bounded in 
C([0,T];L°°). Denote 

W^.= R h V (\a%\ 2 - c) . 

From Corollarvl331 W% G C([0, T];H S+2 ), and VW A £ is bounded in C([0, T];H S 
In particular, Lemma 12.21 shows that Wl is bounded in C([0, T]; L°°). 

From Lemma 12.71 and Arzela-Ascoli's Theorem, for a subsequence h! of 

K 

(3.11) v £ h , -> v e and a|, -» a e in C([0, T]; iff Q ' c ), as /i' -> 0, 

for any s' < s — 1. Moreover, we have v £ ,a £ £ C w ([0,T];X s ). We can then 
pass to the limit in all the terms in ()3.5|) . except possibly the Poisson term, 
that is, the right hand side in the equation for vi . 

To claim that (v £ ,a £ ) solves (|3.1|) - (|3.2|) . we introduce the Poisson poten- 
tial 

VZ-qA^Gn H\ 2 -c). 
Then ()3.5() can be rewritten as: 

( 8 t v £ h + J h {(v eik + v £ h ) ■ VJ h v £ h ) + W pcrt + VV£ = 0, 



(3.12) 



dta £ h + J h ((v eik + v £ h ) ■ VJ h a £ h ) + ^a £ h V ■ v{ = i^AJ%a%, 



AVZ = qG h (\a%\ 2 -c). 

A subsequence of converges in V to some W £ G L°°([0, T] x W 1 ). Since 
V x Wfi = for every h G]0, 1], we deduce that V x W £ = 0. We infer that 
there exists V £ such that W £ = W £ (see e.g. Prop. 1.2.1]), and we note 



(3.13) V 2 V £ G C w ([0,T];H s ). 

On the other hand, Corollary 13.51 and Fatou's lemma imply that |o e | 2 — c G 
L°°([0,T];L 2 ). To prove that (v £ ,a £ ) solves (pTTjl - lpO]! . we now just have 
to check that AV £ — q (\a £ \ 2 — c) =0. We proceed in two steps: first, we 
prove that this quantity is a function of time only. Then, since it is in 
L°°([0, T]; L 2 ), we conclude that it is necessarily zero. We have 

|| V {AV £ - q (K| 2 - c)) || i2 < liminf ||V {AVf t - q (\a £ h \ 2 - c)) \\„ • 
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The last quantity is equal to: 

kl||VJ 1/h (|a|| 2 -c)|| L2 . 
This goes to zero with h, since |a|| 2 — c is uniformly bounded in L^L 2 : 



\VJ 1/h [\al\ 2 -c) 



< 

\l 2 ~ 



V(AK £ -g(|a £ | 2 



< 



I? 



h\\T{\ 



n £ \ 2 
a h\ 



c)\\ L 2<h. 



We infer that 

c))=0, 

that is, AV^ — g (|a £ | 2 — c) is a function of time only. We conclude that 
(v £ ,a e ) solves SH^ -KTh. 

We prove additional regularity for (v £ , a £ ) by showing that {v £ h — v £ , a 6 h —a e ) 
(and not a subsequence) goes to zero in L°°([0, T\; X s+2 x X s+1 ). We will 
use: 

Lemma 3.7. Let ip <E H 1 . Then: 

(1) \\Ji/h¥\\ L 2 -^0 as h^O. 

(2) - J h )p>\\ L , + J 2 )(^|| L2 < 2/ i ||V</>|| i 2. 

(3,) T/iere exisfe C > suc/i i/iaf /or a// /i e]0, 1], ||i?/ l V 2 || L 2^ L 2 ^ C. 

Remark 3.8. Note that in the first point, p> is supposed to be independent of 
h. Otherwise, the conclusion needs not be true, which is easily checked by 
considering tfh{x) = h n / 2 U(hx), where U G S(M n ): || Ji/hfh\\ L 2 = \\j^U\\ L 2, 
is independent of h. In the second point, ip may of course depend on h. 



Proof. For the first point, we write 

||<A/hHL2 = 



L 2 



and we conclude with the Dominated Convergence Theorem. Next, we have: 



lL2 = n(l - 3 (/*)) <h\\(l-j (h^)) |£|£|| L2 , 

since the function 1 — j(h£) is supported in {|£| > l/h}. The second term 
in the second point is treated similarly. The last point follows from the fact 
that the symbol of the Fourier multiplier A^V 2 is bounded. □ 

Denote (w^, d^) := (v £ h — v £ , a 6 h — a £ ), and for s > re/2 + 1, introduce 
Ph(t) ■= KWIIr- + l|V<(i)|| HS + \\d%(t)\\ HB . 
By construction, p e h (0) = 0. Write the equation for {w £ h ,V w £ h , d £ h ) as: 
d t w £ h + {v cik +v £ h ) ■ Vw £ h + w £ h ■ Vv £ = 

= -R h V(I- J 1/h ) (\ a £ h \ 2 - \a £ \ 2 ) 
d t Vw £ h + (v cik +v £ h ) ■ V 2 w £ h + w £ h • V V + Vw% ■ Vv £ h 



(3.14) 



+ Shi 



+Vv £ ■ Vw £ h 



-R h V 2 (I- J x 



/ii> 



,£ |2 



d t d £ h + w £ h ■ Va% + (v eik + v £ ) ■ Vd% + - (d £ h V ■ v £ h + a £ V ■ w 



.£ 

2 



GA<F h + TF h , 
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where the source terms are given by 

St =v\ ■ Vv% - J h (vt • VJ h v%) + q*T x VJ x/h (|a e | 2 - c) 



£ £ h =v\ ■ Val - J h {vl ■ VJ h a%) - i £ -A(I - 4)a%. 



The error term may seem to involve too many quantities (too much 
regularity), compared to the classical approach explained for instance in 
The usual approach would consist in estimating (w|,d|) in I? 
only. We cannot get such estimates because of the Poisson term in S e h : we 
can prove it goes to zero in X s , but not in I? . We proceed in two steps: 

(1) We show that we can apply Gronwall lemma for r e h (t\ with sources 
terms S e h and Yf h . 

(2) We show that these source terms go to zero with h in the norms 
involved at the first step. 

To estimate the first term of r e h , integrate in time the first equation in (|3.14| , 
and use Corollary 13.51 

K(*)IU»</' \\VwUr)\\L°°dT + I \\wUT)\\ L °°dT 
Jo Jo 

+ [ \\R h V(I-J 1/h )(\al\ 2 -\a £ \ 2 )(T)\\ Lca dr+[ ||^(r)||^dr. 
Jo Jo 

Estimate the third term of the right hand side thanks to Lemma 12.21 Corol- 
lary E3] and the last point of Lemma 13.71 

\\R h V(I-J 1/h ){\at\ 2 -\a £ \ 2 ) (r)\\ LOO 

< \\VR h V(I - J 1/h ) (K| 2 - K| 2 ) (r)\\ Ha ^ 

<\\(I-J 1/h )(\ai\ 2 -\a £ \ 2 ) (t)\\ hs 

<||(K| 2 -K| 2 ) (r)\\ H . 
<||(a|-a £ )(r)||^. 

Using Sobolev embedding for the term in Vw|, we end up with: 

(3.15) K(*)[|r~ < fri{r)dr+ f \\S £ h (r)\\ L ^dr. 

Jo Jo 

Now estimate the H s norm of Viuf . From the second equation in (|3.14j) . 

^||A S V<||| 2 =2Re(A s V^,^A s V^) 

< |Re (A s Vwl A s (vi • V 2 ^)> | + |Re (A s W h , A s (w 6 h ■ V 2 v £ ) 
+ Phit? + Pl(t) \\RhV 2 (I - Ji/h)A s {\a%\ 2 - |a £ | 2 ) \\ L2 
+ P £ h (t)\\VSi\\ Hs . 

Write the first term of the right hand side as: 
Re (A s Vw%, A s (vl • V 2 w%)) = Re (A s Vwl v% • V 2 A s wi) 

+ Re (A s Vwl A s (vl ■ V 2 <) - v\ • V 2 A S ^> . 
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Integration by parts and Kato-Ponce estimates ()2.2j) yield: 

|Re(A s V<,A s K-VH))|^^(*) 2 + PU*)l|V^IUoc||V 2 <ll^ 

+ pirnVv^H.-xW^wlW^ <p £ h (t) 2 , 

where we have used Corollary 13.51 and Sobolev embeddings. Similarly, 
Re (A s Vw £ h , A s (w £ h • VV)) = Re (A s Vw £ h , w% ■ V 2 A s v £ ) 

+ Re (A s Vw £ h , A s (w £ h ■ V V) - w £ h ■ V 2 A s v £ ) , 

and: 

| Re (A s Vw £ h ,A s (w £ h • V 2 t> £ )}| < p%[t) 2 + p £ h (t)\\Vw £ h \\ L o°\\V 2 v £ \\ H s-i 

+ P^(i)||v<|| ffs -i||vV|| L oc < P i(t) 2 . 

We also have 

||i4V 2 (/ - J l/h )A s {\a £ h \ 2 - \a £ \ 2 )\\ L2 < \\A° (\a £ h \ 2 - \a £ \ 2 )\\ L2 < p £ h (t), 
and we infer: 

(3-16) ^[iVufllk < p £ h (t) 2 + p £ h (t) \\VS £ h \\ H3 . 

Proceeding similarly for d £ h , we find: 

(3-17) j t \\d £ A<Pi{t? + P £ h{t)\\^ h \\Hs- 

Summing over (|3.15|) and the time integrated Equations ()3.16(l and (|3,17|) . 

we complete the first task of the program announced above: 
(3-18) t 

Pl(t)< f Plir)dr+ r(||^(r)||^ + ||V^(r)||^ + ||S £ ,(r)||^)dr. 
Jo Jo 

Since 5| G H s , Lemma l'2. 21 implies: 

Pl{t)< f Ph(r)dr+ [\\\VS £ h (r)\\ H s + \\j: £ h (r)\\ H s)dT. 
Jo Jo 

It is an easy consequence of Corollary 13.51 and Lemma 13.71 that we have: 

llV^Hiooij-8 + ->• as h ->■ 0. 

We infer from Gronwall lemma that p £ h — > as h — > 0, uniformly on [0, T]. 
Therefore, we have: 

(v £ ,a £ ) G C{[0,T};X s+1 x X s ) ; |a £ | 2 - c G C([0, T}; L 2 ), 

and the existence part of Proposition 13.11 follows by a bootstrap argument 
(to prove the extra smoothness). 

Uniqueness follows from the above computations: up to changing the 
notations, we have the same estimates as above, with now S e = S e = 0. 
Uniqueness then follows from Gronwall lemma. 

To see that there exists <j) £ such that v £ = V(f> £ , apply the curl operator 
to the equation satisfied by v £ (|3.1|) . Energy estimates then show that 
V x v £ = 0. We conclude thanks to Prop. 1.2.1]. 
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Before being more precise about the properties of 4> £ (we already know 
that V<j) £ G C([0, T]; X°°)), we examine the Poisson potential V£. We have 

AV £ = q{\a £ \ 2 -c) E C([0,T}; H°°). 

We infer from Lemma 12,81 that 

Fy* (W p £ ) GC^T]^ 1 ). 

We deduce W £ G C([0, T] x R n ), and Riemann-Lebesgue lemma implies 
that 

VV £ (t,x) 0. 

So far, we have worked with VV^ only, and we know that it is smooth. At 
this stage, V £ is determined up to a function of time only. The condition 
VI(t, 0) = fixes the value of that function, and yields a unique, smooth, 
Poisson potential (so far, only its gradient was unique). As announced in the 
introduction, we explain why we cannot (in general) impose the behavior 

(3.19) V £ (t lX ) — » 0, 

instead of V £ (t,0) = 0. We know from Lemma EH1 that for all t G [0,T], 
W £ (t, ■) G L p (R n ), for p > 2n/(n - 2). We can then apply Lemma l2~Tl only 
when n ^ 5. The following example shows that in space dimension n = 3, 
we may have Vf(x) — > 0, A/ G and f(x) — > +oo: 

>oo \ x \~ *°° 

f(x) = log (l + \x\ 2 ) , x G M 3 . 

Note also that in the case c G L 1 (M n ) discussed below, we have the additional 
property AVI G C([0,T]; L 1 n -ff 00 ), which makes it possible to impose 
(j3.19j) . Back to (p £ , we have: 

v (d t f + v^ cik • v<f + ||V0 e | 2 + y pert + v^j = o. 

We infer: 

d t <\> £ + V0 cik • v<f + i|v<f | 2 + y pcrt + v £ = f, 

where F = F(t) is a function of time only. In the above equation, all the 
terms are uniquely determined, except dt(p £ and F. Imposing 0f t=o = 4>o-, 

and replacing e with (j) £ + J Q * G(t)cLt if necessary, we may assume that 
F = 0. This condition fully determines </> e . This completes the proof of 
Theorem 11.31 

4. Convergence as e — >■ 0: proof of Theorem 11.61 

In this section, we prove Theorem ll.6l First, the existence of (a, 4>) solving 
(|1.8|) follows from the proof of Theorem II .31 since Ijl.SJI is nothing but ()1.7|) 
with e = 0. Denote 



= v £ — v = Vcj) £ — V(j) ; w £ a := a £ — a. 
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The pair {w £ v ,w £ a ) solves a system similar to ()3.14j) : 

d t w% + w £ v -vv + (v cik + v £ )-vw £ v + v (v £ - V p ) = 0. 

d t w e a + w% ■ Va + (v cik + V s ) ■ Vw £ a + ^ «V • v £ + aV • <) = i|Aa e . 
A(V p £ -V p )=g(K| 2 -|a| 2 ). 

v ( f p - v p) -> as M -> 00 ; <L = ° ; w e a \ t=0 = r £ . 

Let s > n/2 + 1. Mimicking the computations made in Section f3,3( we find: 

| (iK(i)ii^ +1 + ikwiii.) < \wmw + 

+e\\Aa £ \\ H s\\w £ a (t)\\ H s + \\V (V^ -V p ) |U*+i|K(t)|| xs+1 . 
From Theorem 11.31 

||Aa e || L§?/ / s < 1. 

To estimate the term corresponding to the Poisson potentials, write: 
II V (V £ - Vp) < || V [V £ - V p ) \\ L ~ + || A {V £ - V p ) \\ H s. 
We estimate the first term thanks to Lemmas 12.11 and 12,31 we have 
A {V £ -V p )=q (K| 2 - |a| 2 ) = q {\a £ \ 2 -c + c- |a| 2 ) . 

Therefore, A (V £ - V p ) is bounded in L^H S . We see from (l3~13j) that 
dj k (V £ — V p ) is bounded in L^H S for every pair (j, k). Lemma l2~T1 (with 
p = 2) shows that there exists a function 7? (t) of time only such that 

Si^-^p) (t,-)-7l(0e^(K n )- 

On the other hand, for all t G [0,T], 

^(v p £ -Vp)(^) ^ 0. 

. . In 

Therefore, 7?(i) = 0, and S,- (V p £ - V p ) (t, •) G L^(M n ). The critical 
Sobolev embedding then shows that 

l|v(v p £ -v p ) \\ l ^ < ||A (v £ - f p ) || L2 . 

Along with Lemma 12.21 fwith g = 2n/(n — 2)), this yields: 

l|v(v p e -v p ) ||l~<||A(V p £ -V p ) 11^, 

and we have: 

l|V (V p e - V p ) || xs+1 < ||A (V p e - V p ) \\ H s < Ww £ a W H s. 

We infer: 

j t [w<mw + w<{t)\? H s) < w<mw + \\<mh+4<®\\H.. 

By assumption, 

IK(0)||a-+i + K(0)||j. = ||r £ ||^ — >0, 

e— >0 

and we conclude with Gronwall lemma: 

IKIU-xs+i + ||u£||z,«>ff« < e+ ||r e ||# s . 
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The strong convergence of the quadratic quantities described in Theorem ll.6l 
follows easily. Note that a similar convergence has been obtained by P. Zhang 
[25] . when V^xt = = ao (hence cj)^ = (3q) and c G L 1 (M ri ). The conver- 
gence in |23j is proved is a weaker sense though (in the sense of measures), 
due to a different technical approach based on the use of Wigner measures. 

To conclude this section, we note that one must not expect ae l ^ +< ^ cik ^ e 
to be a good pointwise approximation of u £ = a £ e l ^ + ^ cik ^ £ . We have: 

yf _ ae i(0cik+</>)/e _ (f e «W>cik+</> E ) / 'e _ ae i(<f>cik+<f>)/£ 

= (a £ - a) e^ik+^/e + ae <W* (j<t> e h _ . 

The first term is 0(e) in L 2 n L°° (we avoid differentiation because of rapid 
oscillations). The modulus of the last term is of order 

\a\ ■ 

Note that our results do not allow us to estimate the argument of the sine 
function. Formally, it should not be smaller than 0(1) in general, so we must 
not expect ae l ^ cik+ ^^ £ to be a good approximation for u £ . To have a good 
approximation, we would have to compute the next term in the asymptotic 
expansion for (a £ ,(p £ ) as e — > 0. We leave out this question at this stage 
here, because we do not have completely satisfactory answers for that issue, 
and resume this discussion when c G L 1 (R n ) below, a case where we have 
more precise information at hand. 

Part 2. Subquadratic eikonal phase 

We now allow the external potential and the initial phase to have qua- 
dratic components. After some geometrical reductions, the analysis boils 
down to the previous one. This reveals some differences though: for in- 
stance, even if | a,g | 2 — c G L 2 , one must not expect \u £ (t)\ 2 — c G L 2 for 
t > 0. 

Assumption 2. Recall that n 3. 

• External potential: V ext G C°°(1R x M n ) writes 

V cxt (t,x) = Vq Uad (t,x) + Vp ert (t,x), 

where V^ uac j G C°°(IR x R n ) is a polynomial of degree at most two in x 
(V 3 F quad = 0), and W pcrt G C(M; H°°). 

• Doping profile: it is a short range perturbation of a constant. For simplic- 
ity, we assume that this constant is 1; 

c = l + c, where c G L 1 n H°° . 

• Initial amplitude: it has the following expansion, 

a%(x) = a (x) + r £ (x), 
where a G X°° is such that \a \ 2 - 1 G L 2 (R n ), and r £ G H°° , with 

\\r £ \\H° — >0, Vs > 0. 

e^O 
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• Initial phase: we have $o G C°°(M ri ) with 

$q(x) = <f> qna >d(x) + 00 (x) , 

where </> qua d *s a polynomial of order at most two, and V^o G H°° . 
Example (External potential). We may take 

n 

V quad (t,x) = ^2^j(t)x 2 j, 
3=1 

an anisotropic harmonic potential with smooth time-dependent coefficients. 
Of course, we may take V peTt G C°°(R; H°°). 

5. The eikonal phase and the associated transport operator 

The generalization of Lemma ll.2l is: 

Lemma 5.1. Under the Assumption^ there exists T* > and a unique 
solution cj) cik G C°°([0,T*] x W 1 ) to: 

(5.1) d t 4>cik + ^|V(/) eik | 2 + Vq Uad (t, X) = ; G ik|i=O = <Aquad • 

This solution is a polynomial of order at most two in x: V 3 i^ c ik = 0. 

Proof. The first part of the lemma was established in 0. Consider the 
Hamiltonian flow associated to 

^| 2 + ^ quad (t,x), 

which yields x(t, y) and y) solving: 

f d t x(t,y) =£(t,y) ; x(0,y)=y, 

( ''" ) \dtt(t,y) = -V s V qfiad (t,x(t,y)) ; £(0, y) = V0 quad (y). 

Following this flow and using a global inversion theorem (see jH] for these 
general results), we construct (fieiki locally in time, but globally in space. 
The idea for the global inversion is to notice that V y x is the identity, plus a 
perturbation which is uniformly bounded in space, and continuous in time 
with initial value equal to zero: there exists T* > such that, for all t G 
[0, T*], y i — ^ x(t,y) is a global diffeomorphism. We denote by y(t,x) its 
inverse. This yields (j) eik G C°°([0,T*] x IP), with 

(5.3) V<f) eik (t,x) = t(t,y(t,x)). 

As a byproduct, the function cf> e i k is sub-quadratic: 9°0 e ik G L°°([0, T*] xR n ) 
as soon as \a\ ^ 2. 

Differentiating ()5.1|) three times with respect to any triplet of space vari- 
ables, we see that ^ = V 3 e ;k solves a system of the form: 

(d t + Vcp cik -V)V =M* ; *| t=0 = 0, 

where A4 G L°°([0,T*] x R n ) (A4 is a linear combination of derivatives of 
order at least two of <^ e ik)- Note that the absence of source term and initial 
datum follows from Assumption [21 Since V^ e ik is given by (|5.3I) . we can 
then use the method of characteristics: setting ty(t, y) = ^(t, x(t, y)) (which 
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makes sense since x(t, •) is a global diffeomorphism) , the above equation 
becomes 

d t ^ = M^> ; $14=0 = ; M G L°°([0,T*] x M n ). 

We conclude with Gronwall lemma that $ = $ = 0. 

Alternatively, one can prove that $ = by an elementary integration by 
parts argument. Namely, since V 2 </> e ik G L°°([0, T*] x M. n ) and $(0, ■) = 6 
L 2 (M n ), we have $ G C([0, T*]; L 2 (R n )) together with the energy identity: 
for all t G [0, T*], 

J \^{t,x)\ 2 dx = 1^ J (A<j> eik (t',x) + 2M(t',x)) \^{t',x)\ 2 dxdt'. 

Hence, again, the desired result follows from Gronwall lemma. □ 

In view of the energy estimates performed in Sectional we will not con- 
sider (|5.1|) . but a nonlinear perturbation of this equation. Indeed, if we try 
to mimic the computations after Lemma l3,4( and after having changed vari- 
ables to work on the characteristics, we have to estimate Dta £ h in I? . From 
the equation, 

\\D t a E h \\ L 2 < (1 + \\v £ h \\ Lao + ||V^|| i2 )(K|| ioo + \\Va e h \\ L2 ) + \\Aa%\\ L2 

+ ll a /(A0eik||L2. 

The last term is new, since now A0 e ik is a non-trivial function (of time 
only). This means that we must not even expect the last term to be finite! 
To overcome this difficulty, we proceed as on the baby model 

d t a + ^aA0 cik = 0, 

where from Lemma 15. 11 A(/> jk is a function of time only. It is convenient to 
introduce the auxiliary function 

a(t,x) = a(t,x)exp Q J A(/) eik (T)dr\ . 

Therefore, it is tempting to replace the condition \u £ \ — 1 G L^?L 2 with a 
condition of the form 



we 2 



Jo A(/) cik (T)dT |2 



1 G LfL 2 . 



Apparently, we have solved the issue mentioned above, but the price to pay 
is that we no longer consider the quantity which is natural in view of the 
Poisson equation. The idea is then to introduce a "ghost Poisson potential" : 



AK £ = ge-/o^e ik (r)dr 



«e2 



/o A0 eik (r)d7 



2 



Ay g = g(Wo A<fek(T)dT -l 
AV~ C = qc. 
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In particular, AV g is a function of time only: V g is quadratic in x, and we 
may choose 

I 1 2 

(5.4) VM, x) = q 1 -^- ( e- ^ A ^(r)<ir _ ^ _ 

Following the idea of [5] , it is consistent to replace F quac j with V^ ua( j + in 
(|5. 1|) . since F g is quadratic and cannot be considered as a perturbation or 
a source term. Even though V g depends on ei k, it is reasonable to try to 
extend Lemma 15.11 Indeed, if we consider the iterative scheme 

2 + WM) = (e-^ 4 &* - l) , 

^eik|t=0 = ^quad, 

with <p^l = 0q U ad , we see that applying Lemma 15.11 inductively shows that 
every iterate is a smooth, sub-quadratic function. We have precisely: 

Proposition 5.2. Under Assumption^ there exists T* > and a unique 
solution c ik G C°°([0,T*] x M n ), polynomial of order at most two in x 
(V 3 cik = 0;, to: 

(5 5) d t ct> cik + ||Woik| 2 + V quad (t,x) = -q^ (e-I^4>M-r)dr _ ^ ^ 

< / ) eik|t=0 = <£quad ■ 

We denote: 

1 /■* 

(5.6) 5 (t) := - / A0 cik (r)(ir. 

z Jo 

Proof. Inspired by Lemma l5.1| we seek directly ^ e ik of the form 
<t>eik{t, x) = t xM(t)x + a(t) ■ x + fi(t), 

where M G M 

n xnO^)i ol £ and /3 G K. Plugging this expression into 
(|5.5|) and identifying the coefficients of the polynomials in x, we find: 

M(t) + 2M(t) 2 + Q(t) = ^"S/o'TrA/Mdr _ ^ ^ . M(Q) = ^ 

a(t) + 2M(t)a(t) + E(t) = ; a(0) = a , 
P(t) + ^\a(t)\ 2 + 1 (t) = ; /3(p) = 0o> 
where 

Vquad(*, a?) = t xQ(t)x + • a; + 7(i) ; <j) qua ,d(x) = t xM Q x + a • x + /Jo- 
Introducing the unknown function = J* * M(r)dr, we see that the equa- 
tion in M can be solved thanks to Cauchy-Lipschitz Theorem applied to 
(M(t),R(t)). Then ct(t) and (3(t) follow by simple integration. □ 

The above proof shows that unless Q(t) = = Mq (a case which boils 
down to Part^), g is a non-trivial function of time. 

The previous result implies that the characteristics associated to the 
transport operator dt + V(j) e Yk ' V present in (|1.7j) can be described very 
easily. 



SEMI-CLASSICAL LIMIT OF SCHRODINGER-POISSON 23 

Corollary 5.3. Letx(t,y) be as defined in (pT2")l . There exist a G C°°([0, T*};'. 
and A G C°°([0, T*]; M nxn (R)), symmetric, such that, for all t G [0,T*], 

x(t,y) = e A ^y + a{t). 

Proof. By Proposition 15.21 V<p e i k (t, x) = M(t)x + a(t) for some symmetric 
matrix M. Since 

d t x(t, y) = V0 cik (t, x(t, y)), 
the result follows by integration. □ 

Remark 5.4. Under Assumption^ V^ c jk is a function of time only, and the 
transport operator c^+V^> e ik-V is trivial. In the above proof, M = 0, and we 
have x(t,y) = y + J^a{r)dT. This relation is reminiscent of Avron-Herbst 
formula (see e.g. 8 ). 

6. Main results 
The analogue of Theorem ll.3l is: 

Theorem 6.1. Let Assumption^ be satisfied. There exists T > indepen- 
dent of e G]0, 1] and a solution u £ G L°°([0,T] x W 1 ) to (fTT Ji -(fO |) . with 

vv p £ (t, x) = v {y; -v s - v~) (t, x ) — ► o, v p e (t, o) = o, 

>oo 

and suc/i i/iai |n £ e 9 | 2 — 1 G L°°([0, T]; L 2 ), where g is given by (|5.6j) . More- 
over, one can write u £ = a £ e l ^ cik+ ^ £ ^ £ , where: 

• a £ G C7°°([0,T] xR")nC([0,T]; X°°), and |a £ e^| 2 -l G C([0, T]; L 2 ). 

• </>eik is given by Proposition \5. l\ 

• <f G C°°([0,21 x R») andV<j> e G C([0, T];X°°). 

• VFe /iai>e i/ie following uniform estimate: for every s > n/2, there 
exists M s independent of e G]0, 1] such that 

\\ a6 \\L°°{0,T;X s ) + ||Ke 5 | 2 - 1^00(0^^2) + \\^4> £ \\l°°{0,T;X s ) < M s- 

£ 

Remark 6.2. We impose conditions on V p , and not on V £ . This is related 
to the arbitrary choice (j5.4[) to integrate the "ghost Poisson equation" (this 
equation introduces additional degrees of freedom), since we will impose 

W H (x) — » ; V~ c {x) — » 0. 

Note that since g is non-trivial, the above result shows that one must not 
expect \u £ (t)\ 2 - 1 G L 2 for t > 0. 

Proceeding like before, we want (a £ ,v £ ) to solve: 

( d t v £ + (v eik + v £ ) ■ Vv £ + v £ ■ Vv cik + Wp Cr t + W p £ = 0, 



(6.1) 



d t a e + (v eik + v £ ) ■ Va £ + ^a E V • (v eik + v £ ) = i^Aa £ , 



AV p £ = g(K| 2 -c), 
together with 

(6.2) VV p \t,x) — ► ; ^ £ (f,0) = ; ^ =0 = V<£ ; 4 =0 = «o- 

With this existence result, we can study the asymptotic behavior as e — > 
of the solution we construct: 
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Theorem 6.3. Under Assumption^ there exists a smooth solution (a, (f>) of 
(ffiTTji withe = 0, such that a,V<f> G C([0,T], X°°), \ae 9 \ 2 — 1 G C([0, T], L 2 ), 
and 

IK - a||i2?fp + ||V(0 £ - 0)||ls?x« — >0, Vs > n/2. 
In particular, the position density and the momentum density converge: 
\u e \ 2 — >\a\ 2 inLZ?H s , and 



elm. ^n e e-^ik/ £ V (u £ e'^^J J — + |a| 2 V0 in L^X S , Vs > n/2. 

Remark 6.4. We slightly altered the usual notion of momentum density, by 
removing first the eikonal phase (f> e ik.- Indeed, we do not prove that 

K| 2 V0 cik — ► |a| 2 V0 cik in LfX s , 

e— >0 

since V^ e ik may grow linearly in x, while |a| 2 morally goes to 1 as |a;| — ► oo, 
hence |a| 2 V0 cik LfX s . 

We first show that the solutions of (|6.1|) exist and are uniformly bounded 
for a time interval independent of e. 

Proposition 6.5. Let Assumption^ be satisfied. Let s > n/2. For all 
M > M Q > 0, there exists < T < T* such that, if for all e G [0, 1], 

(6.3) l|V0 o ||^+2 + ||l«o! 2 " 1|L 2 + IKIIxs+i < M o, 

then the Cauchy problem (|6.1j) has a unique classical solution (v £ ,a £ ) in 
C°°([0,r] x R") such that 

I _e_ai2 



(6-4) IKIlL|?^+a + 



1 



+ IKIIl-a-s+i < M - 



L^L 2 t 



We perform some geometrical reductions so that the proofs of the above 
results follow from Section 

7. Reduction to the first case 

We begin by proving that (|6.1|) is equivalent to a system which does not 
involve the operator v e ik'V, thanks to Corollary 15. 31 Resuming the notations 
of Section |SJ define, for any function / of time and space: 

f(t,y) = f(t,x(t,y)). 

Working with / instead of /, the characteristics associated to v e [ k ■ V are 
straightened so that: 

dj(t,y) = (d t +v eik -V)f(t,x(t,y)). 

The good news for us is the fact that the above change of variable does not 
change the structural properties of (|6.1j) . Indeed, Corollary 15 . 31 implies that 



(7.1) Vf(t,y) = e- A ^Vf(t,y), 

for some symmetric n x n matrix A(t) which is independent of y. 

We are now in position to make precise the fact that the change of vari- 
ables does not change the structural properties. 
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Lemma 7.1. Fix t E [0,T*] and set 5t := e~ A ^V , where A is as in Corol- 
larv \5.c\ The following properties hold: 

(1) For all u E H 2 (R n ) and all v E W^iW 1 ) one has: 

Ke(i5^5tu, u) = ; 2(v ■ 5tu, u) = —((5tv)u, u). 

(2) The Fourier multiplier V(— 51 5t)~ 1 5t is well defined and bounded on 
Sobolev spaces: for all a ^ 0, there exists a constant K a independent of 
t E [0, T*\ such that: 

(3) For all function u : W 1 -> R ; 

u(z) — > 44> u(x(t,y)) — > 0. 

x|^oo I J/| — >-oo 

Proof. By integrating by parts, the first property follows from the fact that St 
is a linear combination of spatial derivatives whose coefficients are constant 
symmetric matrices. The property (2) is immediate using Fourier transform. 
The property (3) is obvious. □ 

Notation. Introduce the operator d by, for all u: [0,T*\ -> <S'(IR n ), 

(du)(t) := e- A(t) Vu(t). 

The difference between the above notation and Lemma 17.11 is that 5t is 
defined for fixed t E [0, T*]. Following what we did in Sectional introduce 

a e := a e e 9 ; V £ = V p + V s + V~ c . 

Since c£ L 1 n H°°, A -1 c is well defined temperate distribution: 

A- 1 c = -^- 1 (|||- 2 JF(c)). 

Setting V peit := V pcTt + gA _1 c, we still have VVjp er t E C(R;H°°), from 
Lemma 12.41 With these notations, Q6.1|) is equivalent to: 

Id t v £ + v £ ■ dv e + v e ■ dv eik + dv pcrt + dv p £ = o, 
d t a £ + v £ • <95 e + ^a £ d ■ v £ = -i^d*da £ , 
d*dV p £ = -qe- 29 (\a £ \ 2 -1). 

Note that the fact that the right hand side of the equation for a £ is skew- 
symmetric remains, from the first point of Lemma l7.ll following the idea of 
E. Grenier |15j . this is crucial in the proof of Theorem ll.3l This is so thanks 
to Corollary 15.31 and would not be if </> e ik was not exactly polynomial. 

Directly from Corollary 15.31 we verify that, for all a ^ 0, there exists Co- 
such that, for all t E [0,T% 

(7.3) C-^K^i,-))!!^ < \\u\\ Hir <C ff \\u(x(t,-))\\ H ., Vu E H a (R n ). 
Similarly, there exists a constant C such that 

C- 1 \\u{x(t, -))\\ LOO < \\u\\ LOO ^ C \\u(x(t, -)\\ LOO , Vn E L°°(lR n ). 

Note that (|7.2|) is very similar to (j3.1|) . The transport operator is simplified, 
but we have two new features: 

• The term v £ ■ <9v e ik i n the equation for v £ . 
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• The factor e 29 in the Poisson equation. 

The latter changes very little computations, since g is a function of time 
only. One can check that the term v £ ■ dv^ does not require a modification 
of the proof given in Section |21 unlike for dtaF , we do not estimate dtv 6 in 
L 2 . Therefore, Theorem Id II follows. Similarly, Theorem 16.. SI follows like in 
SectionHJ up to some slight modifications; like for the proof of Theorem 1(11 1 
replace (a 6 , V £ , V peit ) with (a £ ,V p ,V pert ). 



Part 3. Integrable doping profile 

8. Integrable doping profile 

There are many results concerning the case when the doping profile c is 
decaying at spatial infinity, say c G L 1 (M n ). We refer for instance to [2H] and 
references therein. We restrict our attention to the case n = 3 for simplicity. 

Assumption 3. We consider the case n = 3. 

• External potential: V cxt G C°°(M x M 3 ) writes 

V cxt (t,x) = V quad (t,x) + V pcTt {t,x), 

where V qua ^ G C°°(M x IR 3 ) is a polynomial of degree at most two in x 
(V 3 V quad = 0), and V pcrt G C(R; H°°). 

• Doping profile: c G L 1 (R 3 ) D X°° . 

• Initial amplitude: olq{x) = oq(x) + ea\(x) + erf (x), where aQ,a\,r\ G H°°, 
with 

\\rf\\ H s — >0, Vs > 0. 

e^O 

• Initial phase: we have $0 £ C°°(IR n ) with 

$o0) = <Aquad(^) + 00 

where qU ad * s a polynomial of order at most two, and 4>q £ X 



00 



Our goal is to state a convergence result which is more precise than The- 
orem ESI we shall need some properties of 4> £ , and not only V0 e . This is 
why we change the boundary conditions to solve the Poisson equation: we 
consider 



e 



2 



ied t u e + y Au e = V ext u £ + V £ u £ , (t, x) G M x M 3 , 
(8.1) AV* = q(\u e \ 2 -c), VV £ (t,x) -» and F p e (t,a;) -» as |x| -» 00, 

ujU = aS(x)^W/« 

With these boundary conditions (which are as in [HI ESI f° r instance), we 
can define A -1 as: 

(8-2) A- 1 f = --^ n *f. 

Theorem 8.1. Let n = 3. Under Assumption^ assume furthermore that 
Vp Cr t G C(R;X°°), c G L 1 ^ 3 ), o G L 2 (R 3 ) and (j) G L°°(M 3 ). T/iere exisfe 
< T ^ T* independent of e G]0, 1] and a unique solution u £ G C°°([0,T] x 
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R 3 ) n C{[0,T];H°°) to dSHJ. Moreover, one can write u £ = e e i ^ dk+ * e )/ E , 

• a £ e C([0,T];iP°). 

• <^ e ik «s given by Lemma \5.1\ 

• <f G C([0,r];Z°°). 

• We /iaf e the following uniform estimate: for every s ; i/iere exists 
M s independent of e e]0, 1] swc/i i/iai 

Note that existence and uniqueness for (|8.1j) can be established in a larger 
class of functions, thanks to Strichartz estimates. We refer for instance to 
[2] for the case with no external potential, and simply recall that similar 
Strichartz estimates are available in the presence of a smooth, subquadratic 
external potential ( JOJ) see also [B]). Note also that the term A _1 c can be 
treated as a "nice" linear potential, thanks to Lemma l2.4l and the following: 

Lemma 8.2. The operator A -1 defined by (j8.2j) maps L 1 n L 2 (IR 3 ) to 
JT(L 1 (R 3 )), where T denotes the Fourier transform. Moreover, there ex- 
ists C such that 

\\F (A~V) La < C{\\y\\ L i + \\<p\\v ) , VV> G L 1 n L 2 (M 3 ). 

Uniqueness for H8.ll) follows easily: 

e 2 

ie<9fu £ + yAti E = Kxtn 5 + gA _1 (|u £ | 2 - c) u £ . 

Let u £ and u £ be two solutions in C([0, T £ ]; H°°) of the above equation, with 
the same initial data, for some T £ > 0. Note that the dependence upon e is 
irrelevant, since e > is fixed. The difference w £ = u £ — v £ solves 

£ 2 

ied t w £ + yAiy e = V ext w £ + gA _1 (|u £ | 2 - c) w £ + gA _1 (|u £ | 2 - \v £ \ 2 ) v £ . 
The basic energy estimate yields: 

4lKHi^ll A_1 (KI 2 -KI 2 ) ^LlKIU* 

< IIa-^ki 2 -^! 2 )!!^ |KIMK|| L2 

< (|||u £ | 2 - | V £ | 2 || L1 + ||K| 2 - |t; £ | 2 || L2 ) \\v £ \\ L 2\\w £ \\ L 2 

< (\\u £ \\ L 2 + \\v £ \\ L 2 + \\u £ \\ L ^ + ||v e || L oo) ||«; e [|i2[|u e ||i,2. 

Uniqueness then follows from the Gronwall lemma. 

To prove the existence part of Theorem I8.1| we consider 

d t <j) £ + V0 eik • V<f + \\V<P £ \ 2 + U pcrt + V £ = ; 0f t=o = O . 

(8.3) «9 t a £ + V (<f + cik ) • Va £ + ~a e A (<f + cik ) = i|Aa £ ; a £ 4=0 = a§. 
AU £ = g(|a £ | 2 -c) ; VF £ (t,x) — ► 0andU £ (t,x) — ► 0. 

The geometrical reduction presented in Section [7| makes it possible to trans- 
form the transport operator dt + i> e ik ■ V into dt- Unlike in Section [7[ we 
may keep the term A<^ e i k . Since A0 e j k is a function of time only, and since 
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we work with a £ £ C([0,T]; H s ), the term a e A0 e j k can be treated like a 
perturbative term. 

Since the proof of Theorem 18.11 involves more classical arguments, we 
essentially skip it, so that we can focus our discussion on the semi-classical 
limit e — > 0. 

After the geometrical reduction, (j8.3j) becomes what we would have found 
directly in the case V^ uat j = = </> qU ad> up to terms which can be treated by 
Gronwall lemma. We may for instance resume the approach of Section 
and replace X s with H s . This way, we construct a £ ,v £ £ C([0, T]; H°°). 

To complete the proof of Theorem l8.ll we finally notice that <p £ £ C( [0, T] ; L c 
from Lemma 18.21 and 1)8. 3 j) integrated along the characteristics. 

We can now establish the analogue of Theorem 11.61 with a pointwise 
description. To do so, we introduce the solution to 

5l 



W + V^eik • + -\V(j)\ 2 + Fpert + = ; <j)\ t=0 



(8.4) $ a + V (</> + 0eik) • Va + ^aA (0 + eik ) = ; a| t=0 = a . 
AVp = g(|a| 2 -c) ; VVp(t,x) — ► Oand^^x) — » 0. 

|x|^oo ' |a:|— >oo 

This system has a unique solution (<fr,a) £ C([0,T]; X°° x As pointed 

out at the end of Sectional the triplet (0 e i k > 0, a) does not suffice to describe 
the pointwise limit of u e as e — > 0. This is the reason why in Assumption |HJ 
we want to know a £ up to o(e) instead of o(l) only. Consider the linearized 
system: 

dt<h + V(0 eik + (/)) ■ V0! + V = ; 0!| t=o = 0. 

d t b + V (0 eik + 0) • V6 + i&A (0 cik + </>) + 

(8.5) ! i 

+V0i • Va + -aA0i = -Aa ; b| t=0 = ai. 

AF = 2g Re (ao) ; W(i,x) — ► 0, aadV(*,a;) — ► 0. 

\x\^oo »oo 

It has a unique solution (0 1} 6) £ C([0,T};X°° x il 00 ). 

Theorem 8.3. Under the Assumption^ the solution to Q8.1|) can 6e ap- 
proximated at leading order by ae J ^ >1 e i ^ cik+ ^^ //£ : 



Qg^igiOeik+^/e 



— > as e — ► 0. 

LS?(L 2 nL°°) 



Remark 8.4. In general, 0i is not trivial provided that ai ^ 0, and the 
amplitude of u £ is, at leading order, ae 1 ^ 1 . This phenomenon is due to 
the fact that from the point of view of geometric optics, ()l.lj) - (|l,3j) (or 
Q8.1[) ) is supercritical: to describe the exact solution at leading order as 
in Theorem 18.31 it is necessary to know its initial data up to o(e). This 
phenomenon may lead to instability results as in f3j: modifying Oq at order 
yfe for instance, affects the solution u £ at order 0(1) for times of order yfe. 

Sketch of the proof. The idea is to resume the approach of Section |IJ Set 

wi = V (cjf — 4> — £0i ) ; w £ = a £ — a — eb. 
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Proceeding as in Section^ we find, for s sufficiently large: 



K s+i -r \\w a \\L°°H» ^ t~ 2 
As above, we infer an I? estimate for w%: 

H«^[U°?ff*+i + Wo\\l™h° < e 2 + e\\r\ \\ H », 
and directly from the equation, 

\\<p e -(f)- £^i|Uoo ( [ 0)T ] x]R 3) < e 2 + e\\rl\\ H s = o(e). 
We conclude: 



<e|6| + K| + \a\ 



sin 



6 £ — — £</>i 



<e|&| + |5£| + |a| x (1). 
The result follows by taking the L 2 or the L°° norm in space. □ 
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